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ABSTRACT 

We report the discovery of a new ultraluminous X-ray source (ULX) in the nearby galaxy NGC 4244 
from Chandra archival data. The source, IWCA J1216. 9-1-3743, is one of the least luminous and softest 
ULXs discovered so far. Its X-ray spectrum is the best available for a representative of the soft ULXs, 
a class of sources recently discovered by Chandra and XMM-Newton. IWGA J1216. 9-1-3743 appears 
point-like in the Chandra image and has a complex spectral shape: a multicolor disk model, suitable 
for brighter ULXs, is inadequate for this source. IWGA J1216. 9-1-3743 spectrum is heavily absorbed 
{Nh 1 — 4 X 10^^ cm^^), and very soft. The best-fit power-law model gives F 5 and implies a 



luminosity i(o.5- 



10 koV) 



lO'^-' erg s-^ A comparison with previous detections shows that, despite the 



variability displayed by the source during the Chandra observation, IWGA J1216.9-f 3743 count rate, 
spectral shape and absorption are practically unchanged over a 9-year period. 

We performed also deep optical imaging of the field containing the X-ray source, and found a possible 
R ~ 23.7 counterpart. 



Subject headings: X-rays: galaxies — galaxies: spiral 
individual(lWGA J1216. 9-1-3743 ) 

1. INTRODUCTION 

X-ray observations of nearby galaxies have shown that 
their X-ray emission comes from different classes of 
sources: X-ray binaries, supernova remnants, hot inter- 
stellar medium in addition to possible background active 
galactic nuclei (AGNs; see e.g. Fabbiano 1989 for a re- 
view). Some of the point-like sources appear to radiate 
well in excess of the Eddington limit for a IMq object, 



galaxies: individual (NGC 4244) — stars: 



with inferred luminosities in the range 10'*^ — lO""^ ergs 
s^^ . These objects are often referred to as ultra lu- 
minous X-ray sources (ULXs) or intermediate-luminosity 
X-ray objects (IXOs), inasmuch their luminosity is in be- 
tween those of "normal" X-ray binaries and AGNs. Up to 
now about 90 ULXs have been detected in more than 50 
galaxies (see Colbert & Ptak 2002 for a recent catalog). 
ULXs seem to be preferentially located in the outskirts 
of the host galaxy, although some are found in the cen- 
tral regions (e.g. IXO 95 in NGC 6949). As pointed out 
by Colbert & Ptak, starburst galaxies contain a compar- 
atively large number of ULXs, but spirals seem not to be 
favoured hosts with respect to ellipticals. In a very recent 
paper, moreover, evidence was found for a ULX associated 
with a globular cluster in NGC 4565 (Wu et al. 2002). 
ULXs are often variable sources and their X-ray spectrum 
may be quite soft (e.g. Makishima et al. 2000). No cer- 
tain optical counterpart has been found for these sources 
yet, implying an extreme X-ray-to-optical flux ratio. Very 
recently PakuU & Mironi (2002) reported the presence of 
emission nebulae at the position of some ULXs and sug- 
gested that they could actually be related to the episode 
which led to the ULX formation. 



Thus far no generally accepted model has been pre- 
sented to explain the huge energy output of ULXs. Assum- 
ing that the Eddington limit is not exceeded and emission 
is isotropic implies that ULXs are powered by accretion 
onto a « 50 — 100 Mq black hole. Present evolutionary 
scenarios do not preclude such a possibility and the fact 
that these are rare sources may be compatible with the 
mean number of ULX per galaxy being only 1-2 (Colbert 
& Ptak 2002). Alternatively ULXs may be interpreted 
as conventional black hole binaries {Mbh ^ 10 Mq) with 
modest beaming coming from a coUimated X-ray emission 
{b 0.1, King et al. 2001) or to a relativistically beamed 
emission. In this picture the Galactic analogues of ULXs 
should be the microquasars, such as GRS 1915-1-105 and 
GRO 1655-40. 

This paper presents a multi-wavelength study of a new 
ULX, IWGA J1216.9+3743, found in NGC 4244. This 
source is one of the 16 pecuhar iZ05'^T' PSPC sources se- 
lected for their extremely high X-ray-to-optical flux ratio 
(Cagnoni et al. 2002). IWGA J1216.9+3743 is a bright 
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> 10-13 erg cm-2 s'^) WGACAT (White, 



Giommi & Angelini 1994) source with blank fields, i.e. 
no optical counterparts on the Palomar Observatory Sky 
Survey to 0=21.5. The extreme fx/fv ratio that fol- 
lows is incompatible with all major and common classes 
of extragalactic sources, including normal quasars, AGNs, 
normal galaxies and nearby clusters of galaxies (Maccac- 
aro et al. 1988). Possibilities for the nature of these 
'blanks' (Cagnoni et al. 2002) include: (a) Quasar-2s, i.e. 
high luminosity, high redshift heavily obscured quasars, 
the bright analogs of Seyfert 2s; (b) Low Mass Seyfert- 
2s, AGNs powered by a low-mass obscured black hole (i.e. 
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obscured Narrow Line Sy 1); (c) AGNs with no big blue 
bump, e.g. low radiative efficiency flows; (d) Isolated Neu- 
tron Stars (e.g. Treves et al. 2000); (e) 7-ray burst X-ray 
afterglows or fast variable/ transient sources; (f) failed 
clusters, in which a large overdensity of matter has col- 
lapsed but has not formed galaxies (Tucker, Tananbaum & 
Remillard 1995); (g) high redshift clusters of galaxies and, 
most relevant to this paper, (h) ULXs in nearby galaxies. 

Using X-ray archival data [ROSAT and Chandra) and 
the information obtained from optical and IR follow-ups, 
we present strong evidence that IWGA J1216. 9+3743 is 
indeed a ULX in NGC 4244, a B=-18.4 (Oiling 1996) 
edge-on spiral galaxy (Hubble type Scd) at a distance of 
~ 3.6 Mpc (Fry et al. 1999). 

We present the X-ray, optical and IR observations in § 2 
and discuss the results in § 3. 

Errors in the paper represent 90% confidence levels, unless 
explicitly stated otherwise. 

2. OBSERVATIONS 

2.1. ROSAT 

IWGA J1216. 9+3347 was serendipitously observed by 
ROSAT PSPC in November 1991 and the results of this 
observation are reported in Cagnoni et al. (2002); we sum- 
marize here the main findings. During the 9 ks PSPC 
observation (600179n00) a total of 135 net photons were 
collected between 0.07 and 2.4 keV. We extracted a spec- 
trum and fitted it with both an absorbed power-law and 
an absorbed blackbody model. Despite the scanty statis- 
tics and the large error bars both fits agree on the extreme 
softness of IWGA J1216.9+3743 emission (F = 4.90l^:7^ 
for the power-law and AT ^ 200 eV for the blackbody 
models respectively) and on the presence of absorption a 
factor of 10-30 in excess of the Galactic value in this di- 
rection (NHG',a = 1-69 X 10^0 cm^^; Nh ^ 4.5 and - 1.5 
X 10^^ cm^^ for the power-law and blackbody models, re- 
spectively). The source absorbed flux in the 0.5 — 2.0 keV 
band is ~ 1.2 x 10~^^ erg cm~^ s~^, which corresponds 
to an absorbed luminosity of ~ 2 x 10^® erg s~^ and to 
an unabsorbcd luminosity of ~ 10^^ erg s^^ for both the 
power-law and black body models at NGC 4244 distance. 
ROSAT PSPC lightcurve is consistent with the source be- 
ing constant during the observation (65% probability). 

IWGA J1216. 9+3743 region was observed by the HRI 
instrument on board ROSAT'm June 1996 for ~ 8.5 ks, but 
the source was not detected, consistently with the extrapo- 
lation in the HRI band of the best-fit absorbed power-law. 
The 3(T upper limit on the count rate in a 24" circle is 
0.0023 counts s'^ 

2.2. Chandra 

IWGA J1216. 9+3347 was serendipitously observed 
by Chandra ACIS-S on May 20 2000 during a ~ 

50 ks observation of the nearby galaxy NGC 4244. 
IWGA J1216.9+3743 is 6 kpc off NGC 4244 center and 
is detected point-like source on ACIS-S S2 front il- 
luminated chip with a total of 1948 counts between 0.3 

^ Note that all the points lie within a factor 1.3 
^ E.g. the 42 ks spectrum fitted with an absorbed power law mc 
Nh = 4A6+qH X 10^1 cm-^, r = 4.9llo:25, normalization=3.37lo; 
of freedom (d.o.f.). Errors represent here la confidence. 



and 10 keV with peak position at a = 12'^ 16"" 56.927", 
^ = +37°43'35.89" (J2000). 

We reprocessed Chandra archival data using the CXC 
analysis package CIAO 2.2.1 (Elvis et al., in prep.) and 
CALDB 2.9. The source position determination was im- 
proved by using the new geometry parameter file added 
to CALDB 2.9 and by removing the pixel randomiza- 
tion of ±0.5 pixels introduced by the standard Chandra 
pipehne processing (e.g. Garcia et al. 2001). We cor- 
rected the data for bad pixels and checked for high particle 
background times during the observation. The observa- 
tion is not affected by any significant background flare. 
IWGA J1216. 9+3743 average count rate during this ob- 
servation is (3.92 ± 0.98) x 10^^ photons s~^ between 0.3 
and 10 keV for a total of 1948 counts in the ~ 49 ks of 
net exposure. In a conservative test to confirm whether 
any background flare might affect our results, we have 
excluded the time intervals with background count rate a 
factor 1.2 larger or smaller than the mean value during the 
observation^, as suggested by CIAO 2.2.1 Science Threads 
(http : / /asc . harvard . edu/ ciao/ threads/ acisbackgr ound / ) . 
This would bring to an useful exposure time of 42 ks and 
to 1653 net photons (0.3-10 keV) for IWGA J1216.9+3743 
. A detailed comparison shows that the count rate, spec- 
trum and lightcurves are fully consistent with those ex- 
tracted from the total observation^. 

We extracted a spectrum for the source using a circle 
with radius 17.5" (35 pixels) and a local backgrotmd in 
an annular region with radii of 20" and 50" (40 and 100 
pixels) respectively from which a circle of 15" centered 
on a nearby source was excluded. We binned the spec- 
trum to have a minimum of 30 counts per bin and fit- 
ted it with SHERPA between 0.3 and 10.0 keV with an 
absorbed power-law model and an absorbed disk black- 
body. Both models proved to be a poor representation of 
the data and the results are unchanged when the ACIS 
quantum efficiency degradation is taken into account. We 
tried with more complex models, such as an absorbed 
power-law + disk blackbody, an absorbed power-law + 
a Raymond-Smith thermal plasma and an absorbed bro- 
ken power-law and they arc all acceptable. The spec- 
trum and the residuals to an absorbed power-law model 
are reported in Figure 1 and the results of all the fits 
are reported in Table 1. We also report in Table 1 
the 0.5-10 keV absorbed fiuxes and the 0.5-10 keV ab- 
sorbed and unabsorbed luminosities of the source. It is 
clear from Table 1 that: (i) IWGA J1216.9+3743 spec- 
trum is extremely soft (F ~ 4.901° gg); (ii) simple mod- 
els give an unacceptable fit to the data; (iii) all but the 
broken power-law model predict a high column density 
(~ 1 - 7 x 10^1 cm-2) in IWGA J1216.9+3743 direction; 
(iv) IWGA J1216.9+3743 unabsorbed 0.5-10 keV lumi- 
nosity is > 10^^ erg s~^ (for all but the broken power-law 
model, which could mimic the absorption with a flatter 
low energy slope and for the disk blackbody model, which 
is the statisticcally worst model). 

We extracted lightcurves for the source and background 
between 0.3 and 10 keV using the same circular and annu- 

il gives values perfectly in agreement with those reported in Table 1 
x 10~* photons cm~^ keV~^ and reduced = 1-32 for 39 degrees 
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lar regions used for the spectral analysis. The background 
subtracted lightcurve of IWGA J1216.9+3743 binned 
over 5000 s is shown in Figure 2 and the source displays a 
factor 1.5 variability. A fit with a constant model gives 
a 0.5% probability for constant emission, indicating 
that IWGA J1216. 9+3743 is variable on timescales < 
5000 s. We have confirmed the variability by an inde- 
pendent method using a Bayesian-block analysis to the 
unbinned data which has recently been developed for the 
Chandra Multi- wavelength Project (J. Drake, private com- 
munication; Kim et al. 2002, in preparation). 

2.3. Optical and Infrared 

We obtained deep IR imaging in the K-band (to K ~ 18) 
with NSFCAM at NASA/IRTF on January 2000 and deep 
R-band imaging (to R - 25) with MOSAIC at KPNO 4m 
in February 2001. Wc performed standard imaging reduc- 
tion using IRAF version 2.11.3 and calculated accurate 
astrometric solutions for the data using the digital sky 
survey. 

The Chandra error box (0.67" at 95% confidence level 
for ~ 1000 net photons at 6' off-axis, Kim et al. in prepa- 
ration) contains only one object (at 0.6" from Chandra 
central position), which appears stellar on the R image 
(Figure 3). IWGA J1216. 9-1-3743 counterpart lies in the 
optical extent of NGC 4244, at ~ 6 kpc from the galaxy 
center, and it is thus likely to be associated to the galaxy 
(we exclude a possible AGN in the background in § 3). 

Since the night was not photometric we obtained an es- 
timate of the source magnitude rcscaling from the USNO 
catalog sources in the 17.0 20.0 R-magnitude range falling 
on the CCD. The estimated R-band magnitude is ^ 23.7 ± 
0.5, while the source is not visible in the K-band image 
down to K^ 18. The X-ray over R-band flux ratio, de- 
fined as in Hornschemeir et al. (2001), is ^ 170. 

We searched for radio counterparts from the FIRST and 
the NVSS surveys but none were found within 30" of the 
Chandra position, setting an upper limit of ~ 1 mJy at 1.4 
GHz. 

3. DISCUSSION 

The properties of IWGA J1216.9-I-3743 indicate that 
it is an ULX in NGC 4244. The source absorbed flux 
in the 0.5-2.0 keV band derived from Chandra data (~ 
1.8 X lO^"'^'' ergs cm~^ s^^ ) is consistent within the er- 
rors with that measured with ROSAT (^ 1.2 x 10^^'^ ergs 
cm~^ s~^ ), so the source appears not to have significantly 
varied between the two observations (~ 9 yrs). The unab- 
sorbed luminosity of IWGA J1216.9-I-3743 at NGC 4244 
distance is ^ 2 x 10^^ ergs s""'^ for the absorbed power 
law model Note however, that the 50 ks Chandra obser- 
vation provides evidence for variability on a much shorter 
timescale (^ 5000 s). We also found that the spectral 
shape did not change from the ROSAT observation and, 
although no simple emission model gives a satisfactory 
fit, the X-ray continuum appears to be very soft and can 
be roughly described by a very steep absorbed power-law 
with r ^ 5 but not with thermal emission. A complex 
spectral shape was also found for the stacked spectra ob- 
tained from the ULXs in the Antcnnea galaxies by Zezas 
et al. (2002a) (see their Figure 2): the spectrum presented 
here for IWGA J1216.9+3743 suggests that the complex- 



ity might not be the effect of the superposition of different 
types of ULX spectra, but it is intrinsic to each source. 

Using the X-ray flux at 1 keV from the absorbed power- 
law fit, the R-band magnitude and the upper limit at 
1.4 GHz we computed IWGA J1216. 9+3743 broad band 
spectral indices according to the formula 

log(/(x2)//(a:i)) 
\og(\(x,)l\{x,)) 

We used the canonical 1 keV, V-band (5500 A) and 5 GHz 
points, by extrapolating the R-band flux density to the V- 
band and the 1.4 GHz flux density to 5 GHz assuming a flat 
spectral shape, and we obtain ano > 0.60, aox = 0.23 
and aRx > 0.47. Such values place IWGA J1216.9+3743 
out of the region occupied by AGNs in the auo^ oiox plane 
(e.g. Caccianiga et al. 1999) since not even the most ex- 
treme BL Lac objects can reach such values (e.g. Costa- 
mante & Ghisellini 2002). Taking absorption into account 
would make the situation even more extreme (e.g. an ab- 
sorbing column of ~ 4 X 10^^ cm~^ with Galactic dust 
to gas ratio would enhance the X-ray flux by a factor of 
^ 10 and the optical flux of a factor ^ 6). We can thus 
exclude the possibility of IWGA J1216.9+3743 being a 
background AGN. 

We can also exclude the possibility of 
IWGA J1216. 9+3743 being a supernova remnant for 
the short timescale variability displayed in the Chandra 
observation and by the lack of a significant fading in the 9 
years between the ROSAT and the Chandra observations. 

Even if the best fit model for IWGA J1216.9+3743, a 
convex absorbed broken power-law, is usually used to de- 
scribe the synchrotron peak observed in beamed sources 
such as blazars, the sharp change in the power-law shape 
of IWGA J1216.9+3743 (Fi = 0.22, = 3.85) is diSicult 
to reconcile with the smoothly curving blazar synchrotron 
peak (e.g. Fi = 2.1, Fz = 2.8 measured for Mrk 421 by 
Guainazzi et al. 1999). 

IWGA J1216.9+3743 luminosity falls at the lower end 
of the ULX range. Some emission models implies Lx < 
10^^ ergs (see Table 1). In particular, the absorbed 
power-law plus Raymond-Smith model, and the absorbed 
broken power-law models, which have the lowest x^, give 
a luminosity of only few 10^* ergs s~^ . This would 
make IWGA J1216.9+3743 an ordinary X-ray binary in 
NGC 4244. However, optically thin brcmsstrahhmg emis- 
sion is typical of extended sources, like clusters of galaxies, 
and we found no evidence for a diffuse nature of this source 
in Chandra data and the broken power-law model does 
not have a straightforward physical interpretation. The 
relatively low luminosity of IWGA J1216. 9+3743 may 
still suggest that this source is an X-ray binary with a 
M « lOM© black hole, of the Cyg X-1 type. The X-ray 
spectrum however strongly argues against this possibility, 
being much softer of those of Galactic black hole candi- 
dates (BHCs). BHCs spectra in the high state are ther- 
mal and soft, but peak around a few keV and extend up 
to ~ 10 keV. 

The extreme softness of the X-ray spectrum makes this 
source rather peculiar among the ULXs. In fact, while very 
soft spectra have been already detected from other ULXs 
with ROSAT, successive observations over a wider energy 
range showed a harder component. An example of this 



4 



Cagnoni et al. 



behaviour is MS 0317.7—6647 which has been first associ- acknowledges a CNAA fellowship and A.C. the MIUR. 
ated with an isolated neutron star by Stocke et al. (1995) 
on the basis of its soft, thermal spectrum. The source was 
later identified with a ULX by Makishima et al. (2000) 
when ASCA data convincingly showed a hard tail (F ^ 2), 
possibly associated with a multicolor disk blackbody. Be- 
fore the advent of Chandra the ULXs seen by ROS AT a,nd 
ASCA were consistent with emission from a multicolor disk 
blackbody (e.g. Makishima et al. 2000); Chandra is now 
discovering a new class of steep ULX (e.g. out of the 30 
ULX of Zezas et al. 2002b, 9 have a steep spectrum). 
IWGA J1216.9-I-3743 is the source for which the best 
Chandra spectrum is available so far. In the study of the 
Antennae galaxies Zezas et al. (2002a) find steeper spectra 
for the least huninous ULX: IWGA J1216. 9+3743 seem to 
fit in this picture, interpreted with the possible presence of 
undetected diffuse hot interstellar medium in the proxim- 
ity of the source. Alternatively one could consider that the 
indication of a luminosity dependent spectral shape is due 
to difi^erent accretion properties. The steep ULXs could 
be the "galactic" analog of Narrow Line Seyfert galaxies 
(a steep version of the Seyfert galaxies) which are thought 
to be powered by a lower mass black hole accreting at a 
higher rate compared to "normal" Seyferts. 

Follow-up optical observations allowed us to discover a 
possible counterpart of IWGA J1216.9-h3743 at 0.6" from 
Chandra position. The association is based only on the po- 
sitional coincidence. The Chandra error box is quite small 
(95% confidence radius < 1", Kim et al. in preparation) 
and the background field is not very crowded (see Fig. 3), 
so the possibility of a chance alignment appears unlikely, 
albeit it can not be ruled out on the basis of present data. 
The optical counterpart has i? ~ 23.7; if this source is a 
star in NGC 4244, it has to be a red star or its bolometric 
correction would imply a luminosity not compatible with 
even the most himinous stars known. The most likely pos- 
sibility is that the suggested counterpart is a late type 
giant /supergiant, like a M5 II-III; for this spectral range 
the source luminosity, corrected for extinction, would be 
L » 7 X 10^7 ergs s'^ , about 10* L©. The R-K< 5.7 
derived from the IR non-detection is not stringent and it 
is satisfied by the red stars. 

ULXs are one of the most interesting class of sources 
which are now being investigated by the X-ray satel- 
lites Chandra and XMM-Newton. We presented in this 
paper the best spectrum (^ 2000 photons) available 
so far for a steep-spectrum ULX. Steep ULXs, like 
IWGA J1216.9-I-3743, , appear to be related to the least 
luminous objects of the ULX class (e.g. Zezas et al. 2002a) 
and were essentially unknown before the launch of Chan- 
dra and XMM-Newton. Future investigations are needed 
to understand the emission mechanisms powering these 
sources and to confirm and explain a possible luminosity 
dependence of the spectral shape. 

This work made use of the Digitized Sky Surveys pro- 
duced at the Space Telescope Science Institute under U.S. 
Government grant NAG W-2166 and of the US Naval Ob- 
servatory catalog USNO-Al.O manteined by the European 
Southern Observatory. I.C. thanks Luigi Foschini and An- 
dreas Zezas for the useful scientific discussion. We are 
also grateful to the staffs of NASA-IRTF and FLWO. I.C. 
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Fig. 1. — Spectrum and residuals to an absorbed power-law fit. The residuals show that a more complex model is needed to model 
IWGA J1216.9+3743 emission. 
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Fig. 2. — Chandra 0.3-10.0 keV lightcurve of IWGA J1216. 9+3743 binned over 5000 s. The dashed line represents the mean count rate 
during the observation (i.e. 3.86 x 10~^ counts s~^) 




Fig. 3.— ROSAT PSPC (r = 39") and Chandra (enlarged to r = 2" for clarity) error circles on IWGA J1216.9+3743 positions. North is 
up and east is left in the image. 



Table 1 

Fit to the 0.3-10 keV Chandra spectrum with SHERPA 







Model 




Parameters values 
{Eb, kT in keV) 


xl/d.oI. 


Abs. 
Fx 


Abs. 
Lx 


Unabs. 
" Lx ^ 


Abs 


X 


power-law 


A 1 Q+0-75 
*-iy_0.67 


r = A.90toil K = {3.38toil) X 10-4 


1.38/47 


2.31 


3.6 


21.6 


Abs 


X 


disk blackbody 


n C _|_g 

U.85_o:47 


kT = 0.300Io;o27 K i-l^zfl^ 


1.91/47 


2.19 


3.4 


4.6 


Abs 


X 


(power-law -|- disk blackbody) 




r = 6.88tJ;?8 /fp,„ = (6.9lt|^2) ^ ^(3-4 


1.28/45 


2.11 


3.3 


104.5 


Abs 


X 


broken power-law 


0.169'* 


ri = o.22tEJf3r2 = 3.85tEJ-?^ 

= 0.94 ± 0.04 K = (1.71^0:22) X 10"'' 


0.89/46 


2.15 


3.3 


3.5 


Abs 


X 


(power-law -|-RS) 




r = 4.26t° ifpow = IMt^oil X 10-4 
kT = 0.8351°;°^^ ^bun = O-S^Rrs = (1.53l° j^) x 10-4 


1.05/45 


2.14 


3.3 


9.4 



'^in units of 10^^ cm-^ 

''Between 0.5 and 10 keV in units of 10^^ erg cm-^ 
■^Between 0.5 and 10 keV in units of 10^^ erg s-^ 

^The absorption was fixed to the Galactic value in IWGA J1216. 9-1-3743 direction; if left as a free parameter it tends to zero 
^Similar values are obtained for Raymond-Smith abundance equal to the solar value 



